Classical thermodynamics of gravitational collapse 
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We study numerically gravitational collapse of a spherically symmetric instanton 
particle in five dimensions. We show that the late stages of the process are character- 
ized by a nearly constant "free energy" , the value of which matches (within numerical 
uncertainties) the value obtained from standard black- hole thermodynamics. This 
suggests a purely classical interpretation of the free energy of a black hole. 



INTRODUCTION 



Brane-world scenarios [H, 3 provide a new motivation for studying energy relations as- 
sociated with gravitational collapse. In these scenarios, instanton transitions 0, 0| on the 
brane can be viewed as transport of the topological charge through the brane, much like 
a change of magnetic flux through a superconducting ring can be viewed as transport of 
vortices through the ring. In this picture, a 6'-vacuum [5|,l6[ corresponds to a steady flow 
^! of the topological charge into extra dimensions, and preventing such a flow amounts to a 
r-| I solution to the strong CP problem. The simplest mechanism that can achieve that is an 
energy cost 0, [1] . Indeed, while in the 4-dimensional chromodynamics instantons connect 
■ states that are exactly degenerate in energy (a consequence of the symmetry with respect 
^ ! to the "large" gauge transformations [s], [J), in the presence of extra dimensions that need 



not be so [9l 



■ In one version of the scenario, energy accumulates in an extra-dimensional black hole 

Tol. However, it is not clear from the outset what is the appropriate definition of energy 



^ , in this case. Should we be talking about the usual mass of a black hole, as defined by the 
00 I spatial asymptotics of the metric, or is there a different "free energy" associated with late 
O ■ stages of gravitational collapse? Evolution of a quantum state on the brane is determined by 
exp[i/(t)], where I is the action of matter and gravity (including those in the bulk). Based 
^ ! on this, in the present work, we explore properties of the following quantity 

m-'-f^ (X) 

i.e., minus the Lagrangian computed on the collapsing solution. In a brane-world scenario, 
there is a preferred time coordinate — the proper time of an observer on the brane, and it is 
this time that we use in the definition ([T]). 

We consider evolution of an instanton "particle" in the 5-dimensional (5d) SU{2) Yang- 
Mills (YM) theory coupled to gravity in asymptotically flat spacetime. Although this asymp- 
totically fiat case is not a realistic braneworld (gravity and the gauge field are not localized 
on the brane), the question about time-dependence of F can still be asked; by a "brane" in 
this case we mean a sphere at large radius. Note that the derivative in ([1]) is a full derivative, 
rather than a partial one taken at fixed endpoint values of the canonical coordinates; the 
partial derivative would give us the mass E of the black hole. 

We refer to F as the free energy (of a classical black hole), for the following reason. 
Suppose at late times F{t) becomes a constant, i.e., the action becomes a linear function of 



2 



time. Such a function can be analytically continued to the Euclidean time r = —it, so that 
the Euclidean action is Ie = Ft. Since we are dealing with just a single classical trajectory, 
the partition sum reduces to exp(— F/T) (where T is the Hawking temperature Ty), which 
means that F is the free energy. 

Interpretation of the Euclidean Lagrangian as a thermodynamic potential is known in 
the literature and, as we will see, our result agrees with the corresponding thermal 
expression, Eq. ([3]) below. However, this agreement is by no means obvious a priori: the 
thermal result is obtained by considering the stationary vacuum exterior alone, while our 
result corresponds to a time-dependent metric of collapsing matter. Moreover, during col- 
lapse, F{t) is expected to change, and a priori it is not even clear if it will approach any 
constant value at late times. 

The latter point, however, can be checked directly in the case of spherical symmetry, 
when collapse can be simulated numerically. In this paper, we present results of such a 
simulation. 

Gravitational instability of an initially static instanton particle in the 5d pure YM theory 
follows directly from the classical scale invariance of the corresponding 4d theory and had 
been discussed previously [l3|. Here, we discuss the process of gravitational collapse of such 
objects (we have considered both static and nonstatic initial configurations). We observe 
formation of a black hole (of some energy E) and find that soon after that F reaches a 
nearly constant value F ^ ^E. 

We interpret the factor | as follows. For a 5d Schwarzschild black hole, the Hawking tem- 
perature, computed by the method of Ref. [l2|], is T = h/ATrrQ, where rg is the gravitational 
radius: 

rl = ^G,E , (2) 
and Gs is 5d Newton's constant (our choice of the radial coordinate r will become apparent 



shortly). The entropy Uj, llJ] is 5" = ATi'^rQ/hG^. Then, by standard thermodynamics, the 
free energy is 

-Pthermo = E — ST = -E . (3) 

Note that, while T and S each contain a power of the Planck constant h, the product ST 
does not and so has a finite classical limit. We see that the value ([3]) agrees with the value 
of F reached in the course of our classical simulation. 

For the Schwarzschild exterior alone, Eq. ([T]) would give F = E. We find that the 
difference F — E comes from a thin shell near the horizon. It may seem, then, that F being 
nearly constant at large times is simply a consequence of the slowing down (redshift) of the 
near-horizon dynamics as seen by a distant observer. If the near-horizon evolution were to 
stop entirely, however, we would again have F = E. (This relation holds for any static 
configuration, vacuum or not). This not being the case means that the fields continue to 
change, and their time dependence significantly contributes to F. The situation is similar to 
ordinary statistical equilibrium: macroscopically nothing changes, but there is still motion 
on the microscopic scale. 

In the remainder of the paper, we describe various technical details of our approach. We 
conclude with some comments on the quantum behavior of the system. 
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II. EQUATIONS OF MOTION IN SPHERICAL SYMMETRY 

A. Spherically-symmetric ansatz 

We use the isotropic coordinates, in which a spherically-symmetric metric takes the form 

ds^ = -N^{t,r)dt^ + ^\t,r){dr^ + r^dnl) . (4) 

These are convenient for numerical studies of gravitational collapse because of the absence 
of singularities in the metric functions: formation of a black hole is reflected in vanishing 
of the lapse function N at some point, while the conformal factor \l/ remains finite every- 



where. This observation was used previously in Ref. [15| for studying the evolution of metric 
perturbations after inflation. 

The 5d Schwarzschild geometry in the isotropic coordinates has the form 

'^0 1 / 0\ 2 



ds' = - ^^ i^ dt' + f 1 + ^ ) {dr' + r'dnl) , (5) 



1 + 1 



where vq is the gravitational radius. Note that \E' does not stay finite in the "white-hole" 
region of this geometry — there is a singularity at r = 0. We find, however, that, while the 
exterior metric of black holes that we obtain numerically well matches the r > vq region of 
the metric the interiors are quite different, and in our case \l/ stays finite. 

When the metric is asymptotically fiat, the total energy (ADM mass [l^) of the spacetime 

is 

Etot = -^R' dr^ilU^ , (6) 

where R is some large radius. In the limit R oo, i?tot will be conserved. One should 
keep in mind, though, that, for a general initial state, only part of the mass will collapse to 
form a black hole, while the rest will be radiated away in the form of matter waves. The 
expression ([6]) is the total and does not, in general, coincide with the mass E of the black 
hole. 

Related to the isotropic coordinates are Cartesian coordinates x*, i = 1, . . . , 4, such that 
dr"^ _j_ 7-2^^2 _ (^g^3;«^2_ these Cartesian coordinates, our ansatz for the SU{2) gauge field 
is 

where /i = 0, . . . , 4; a = 1, 2, 3; nj is the unit radial vector, rjlj are the t'Hooft symbols 
and f{t, r) is the dynamical variable. Note that = at all times, which is consistent 
with the equations of motion. 

For a gauge field of the form ([7j), the topological charge 

Q = -7^ f d'xe^^'''F%F\, (8) 



647r2 / 

is the field strength) becomes 

3 /■ „, , I r d 



Q = --Jfif- 2)f'dr = --Jl{f- V')dr . (9) 
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So, any configuration that interpolates between / = at r = and / = 2 at large r has 
unit topological charge. An example is the BPST instanton [3|, 

2r2 

foir) = , (10) 

where A is an arbitrary scale parameter. 



B. Equations of motion 

The part of the Einstein- YM equations that are evolution equations can be written in 
the Hamiltonian form 



L = JL 
N ^2' 



(11) 



V f N' dV 



(13) 



4 ' 

Here = SttG^, dots and primes denote temporal and spatial derivatives, respectively, 

(15) 

and 

V can be thought of as a potential for the gauge field, and S as the field's energy density. 
The instanton (ITOl) is then a "kink" , which interpolates between two degenerate minima of 
V^, / = and / = 2. 

We see that K (the trace of the extrinsic curvature) and p can be regarded as conjugate 
momenta of \1/ and /. The lapse function A^, on the other hand, has zero conjugate momen- 
tum, i.e., it is a Lagrange multiplier. The Einstein equation for it is an ordinary differential 
equation that has to be solved at each moment of time 



The remaining Einstein equations are constraints: the energy constraint 

.V^^ 3 



where is the fiat-space Laplacian, and the momentum constraint 

U,. if 



= T^^^/' • (19) 
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If these are satisfied at the initial time shce, the evolution equations together with Eq. ([T 
ensure that they are satisfied at all times. 

Eqs. ( IT2!) and (IT^ require boundary conditions at r = _R. We choose 

K'{t, R) = p'{t, R) = 0. (20) 

Although this does not automatically enforce the momentum constraint (fT9l) at r = R, for 
a large enough R the constraint holds there to a sufficient accuracy. 

Regularity of the gauge field at r = implies that f{t, r) = O(r^) and p{t, r) = O(r^) at 
all times. Using these in the momentum constraint fll9l) . we find 

X'(t,0) = /(t,0) = 0. (21) 

Together with Eq. ( |20l) . this provides enough boundary conditions for the evolution equa- 
tions. 

We still need boundary conditions for Eq. (fT7|) . The metric ansatz (jlj) leaves us the free- 
dom of a global (r-independent) redefinition of time. We use this freedom to set N(t, R) = 1. 
If the boundary is far away from any matter, at large r is close to its Newtonian limit, 
for which 

N-(t,R) = -?^. (22) 

For large enough R, the difference of this from zero is inessential, and the data presented 
here have been obtained with 

N'{t, R)=0. (23) 
Finally, we list a convenient integrated form of the energy constraint 



Etot = 27r2 [ (£- —K^ ] mh^dr . (24) 



R / Q 

V 

The constraints are not involved in the evolution algorithm but are used to monitor the 
accuracy of the simulation. 



C. The action 

The original action of the model is 

where Ib is the boundary term chosen so as to remove the second derivatives of the metric 



functions [17[. Just as we did for the energy, we need to distinguish between the total action 
(12^ and the part of it, which we call J, that corresponds to the black hole alone. This is 
possible at sufficiently late times, when the black hole and the outgoing radiation become 
well separated. 

When computed on our spherically-symmetric ansatz, the total action becomes 



I tot 



R 

dt / drL , (26) 
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where 



-iVr 6/2-6^ + 12^ . (27) 



A static configuration is one for wliicli \I/ and / at a given instant vanisli. For any sucli 
configuration, one can integrate by parts the term with N' in fl26|) and use the energy 
constraint (1181) to bring the Lagrangian to the form 



ditot _ 
Hr ~ AG, 



R-^ N^^'l^^ . (28) 



In the asymptotically flat case, and \1/ approach unity at large r, and the right hand side 
of Eq. (!28|) becomes minus the energy IQ. This is an expected result. Indeed, for static 
configurations the total and partial derivatives of the action coincide, so the Lagrangian 
equals minus the Hamiltonian. 



D. Initial state 



In general, an initial state is an arbitrary solution of the constraints (fTSll and (fT9l) . 
One class of such states are static states, for which the momenta vanish identically, so the 
momentum constraint is trivially satisfied. To solve the energy constraint for such states, 
we pick a function /(r) for which gravity is weak and solve the constraint equation 



- ^V'^ = ^ + Vj (29) 

by a formal expansion in the parameter ( = k'^/qym- The actual expansion parameter is 
C/X^, where A is the spatial size of the configuration. For example, if /(r) is taken to be the 
instanton (fTOj) . the first three terms in the expansion are 



2C(2A2 + r2) 2C2(27A'5 + ASX^r^ + 32X^r^ + 8 



,6^ 



^ 3(A2 + r2)2 45^2(^ + ^2)4 +••• • (30) 

Retaining a sufficient number of terms allows us to solve the constraint to any required 
accuracy. 

Using a large static instanton as the initial condition offers certain advantages for nu- 
merical work. These are related to the fact that such an instanton is nearly a classical 
solution (it is destabilized only by gravity, which for large A is weak). We find that in this 
case collapse generates relatively little of an outgoing wave, and the total action Jtot nearly 
coincides with the action of the black hole, / ~ /tot- This simplifies the measurement of /. 
[In general, / has to be obtained by integrating over r in fl26l) not from to i? but from 
to some intermediate radius, so that the outgoing wave is left out.] 

Initial stages of the collapse of an instanton with a large A can be studied analytically, 
since at large A the system is in the Newtonian limit. The equations of motion f[TT]) - f[T2|) in 
this limit combine into 

f N' dV 

f-r^^ + wf'-jf- '''' 
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We search for a solution in the form 

^) = r^ + X^{t) + + • ^^^^ 

This describes the slow motion of the instanton size ("modulus") A, together with small 
deformations of the instanton shape that this motion causes. 

Now, in terms of the parameter ( introduced above, A is 0(a/C), while A is 0{(). So, in 
fl3T|) we need N' /N only to order (. Using Eqs. f|T71) and fl30|) . we find 

^ = -^ + 0(C=) = A(^__ + __j+0(C^). (33) 

Substituting Eqs. fl32|) and fl33|) into Eq. fl3T]) and projecting onto the instanton zero mode 
t/'o = dfo/dX, we obtain an equation for A(t), 

The solution that has a maximum size Aq at time t = is 



This allows one, for instance, to determine the time it takes the instanton to collapse (see 
below) . 

We have also considered nonstatic initial conditions with / of the form (|TOl) and the 
momentum p of the form 

where f is a parameter. (This does not quite satisfy the boundary condition fl20l) . but for 
large R the mismatch is small.) The results are similar to those presented below for the case 
v = 0. 



III. NUMERICAL RESULTS 



Following Ref. [15|, to improve stability of the numerical method we add to the right- 
hand side of (HM a term proportional to the energy constraint. (This brings in the flat-space 
Laplacian of \1/.) We choose the unit of length so that the gravitational radius ([2]) is tq = 1 
for E = -Einst = S^^/S'yM i^^^ instanton mass in the absence of gravity). This is equivalent 
to setting ( = K'^/gyM ~ f ^^e only parameter that appears in the equations of motion. 
The energy and the Lagrangian are made dimensionless by measuring them in units of -Emst- 
In these units, Eq. ([2]) becomes 

E = ry (37) 

A change of the spatial coordinate 

r ^ x = (38) 
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FIG. 2: ^{t,r) as a function of r at different moments of time. 



allows us to discretize the equations on a uniform grid (in x), while maintaining a higher 
density of grid points (in r) in the region of main interest r ~ 1. Results are presented, 
however, in terms of the original coordinate r. 

These results are for the static initial state (fTOj) with A = 5. Figures [TH3] show the time 
evolution of the profiles of the three dynamical functions, /, and N. Formation of a black 
hole is signaled by N{r, t) (as a function of r) crossing zero. In the present case, this happens 
at t ~ 28. Note that the rapid evolution shown in the figures follows a prolonged stage during 
which little happens. The duration of that stage depends on the initial instanton size and, 
for large A, can be estimated using Newtonian gravity: using C = I in Eq. ([35]), we find 

t = ^X\ (39) 
which is well born out by our numerical results. 
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FIG. 3: N(t, r) as a function of r at different moments of time. 

At late stages, the metric functions outside the horizon approach the Schwarzschild form 
([5]), while inside the horizon they approach a curious collapsing universe with nearly space- 
independent and \E'. 
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FIG. 4: Dependence of various quantities on time, ditot/dt is the sum of the gravity {dlgnw / dt) 
and matter (dlmat/dt) contributions, each of which is the integral over r of the corresponding term 
in Eq. (I27p . "Horizon^" is the square of the radius at which N{t,r) crosses zero. 

Various time-dependent quantities are plotted in Fig. HI "Energy Constraint" refers 
to the right-hand side of the constraint equation (^^, and its deviation from i^tot is a 
measure of the accuracy of the simulation. We see that the constraint deteriorates around 
t = 40 (as expected from the presence of large gradients), but before that the free energy 
-^tot = —dItot/dt reaches a nearly time-independent value of about 0.3. The energy E of the 
black hole at late times can be computed from the horizon radius using Eq. (137|) and, as 
seen in the figure, at late times is about 0.9. 

Fig. [5] shows where the free energy is accumulated. The growth from —0.6 to 0.3 in the 
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FIG. 5: Accumulation plots for the action (j26|) at different moments of time. (These show the 
negative of the Lagrangian contained in the ball of radius r, as a function of r.) 

region outside the horizon is the contribution to Ftot from the static Schwarzschild exterior 
and, as expected, closely matches the energy of the black hole. Note that while an outgoing 
wave with a crest at r ~ 10 is clearly seen in Fig. [H it practically does not contribute 
to Ftot- In other words, in this case, Ftot and F, the free energy of the black hole alone, 
nearly coincide. Combining the results of Figs. [Hand El we find that at late times F ^ ^E. 
Moreover, we see that the difference between E and F comes from a thin shell near the 
horizon (the drop at r ^ 1 in Fig. [5]). 
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IV. CONCLUSION 

The counting of powers of h, described in the Introduction, shows that, unlike the tem- 
perature or entropy, the free energy of a black hole remains finite in the classical limit and 
so may have a purely classical interpretation. Our results indicate that this interpretation 
is furnished by Eq. ([1]). At present, our evidence for this is purely numerical and confined 
to a specific type of black hole in a specific 5d theory. However, if the correspondence turns 
out to be generic, it would be interesting to see if it can be reproduced by an analytical 
argument. 

We conclude with a remark on the quantum theory of the system described by the action 
fl2^ . In spherical symmetry, quantization of gravity is straightforward: gravity has no 
independent degrees of freedom and simply follows quantized matter. As we have seen, 
in the Newtonian limit, gravitational collapse of an instanton is described by the ordinary 
differential equation for the modulus A. (The Newtonian limit applies for A ^ tq.) The 
corresponding Hamiltonian is that of a "fall to the center" : 

n2 ACM , , 

where 11 = MA is the momentum conjugate to A, and M = 2Ei^st = ^QT^^/gyM- The overall 
normalization of the Hamiltonian can be restored by noting that, while the first nontrivial 
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term in the expansion (130|) gives E'tot = -E-inst, the second gives a correction, which should 
equal the potential term in (HOl) . 

As well known, a quantum problem described at large A by Eq. (l40l) is in one of the two 
different regimes, depending on the value of the coupling 

o='SMl. (41) 

For 0" > |, there is an infinite number of bound states with small binding energies. For 
0" ^ 1, these states can form superpositions (wave packets) that move nearly classically. 
From the relation (HTl) . we see that this classical limit corresponds to large values of -Emst, 
while the critical cr = ^ corresponds to -Einst being of the order of the 5d Planck mass. 

This work was supported in part by the U.S. Department of Energy through Grant DE- 
FG02-91ER40681 (Task B). 
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